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Abstract 



I give an overview of some of the most used measures of entanglement. To make the 
presentation self-contained, a number of concepts from quantum information theory are 
first explained. Then the structure of bipartite entanglement is studied qualitatively, before 
a number of bipartite entanglement measures are described, both for pure and mixed states. 
Results from the study of multipartite systems and continuous variable systems are briefly 
discussed. 
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Preface 



"What's the name again of this quantum effect when one particle on one side of the universe 
magically affects a particle on the other side of the universe?" , I asked while walking out 
the door to the corridor. I was on my way to a group of teenagers who were waiting to hear 
how it is to study physics at the university. In case someone asked why I found physics 
so exciting I could always mention this quantum, non-locality, spooky action-at-a-distance 
thing. If I only remembered the name of it. . . "That'd be entanglemenf\ a graduate student 
next-door suggested, using the English term. Yes, that was the word lost in the back of my 
mind. I had actually never heard a Norwegian name for it. 

This happened one of the last days of October 2003. About a month later it was 
decided that my master's thesis would be about this very action-at-a-distancc thing, or 
more precisely, about how to quantify entanglement. From hardly remembering the name 
of the game to understanding how an appropriate measure needs to behave has been a long 
and interesting journey. In writing this text I have tried to keep in mind my own state of 
knowledge as of October 2003. With that I hope that any 5th year physics student will 
find it sufficiently understandable and self-contained as to be able to grasp the main points. 
Whether I have succeeded or not is for others to decide. 

This thesis is the final work leading up to the degree of Master of Science in Applied 
Physics and Mathematics at the Norwegian University of Science and Technology (NTNU). 
It was written during the 2004 spring term. 

I would like to thank Jan Myrheim and Stein Olav Skr0vseth for taking the extra effort 
to find a topic in quantum information theory for my thesis, and for supervising it. I would 
also like to thank Viktor H. Havik and Henrik Tollefsen for interesting discussions about 
entanglement. 
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Chapter 1 

Introduction 



Quantum entanglement was first viewed as a curiosity when it was pointed out in the 
attempt by Einstein, Podolsky and Rosen |EPE,35j to show that quantum mechanics could 
not be a complete theory. When John Bell many years later published his now well-known 
theorem jBel64l , it was hardly noticed in the scientific community. Bell had showed that for 
a theory based on local hidden variables, certain correlations were upper bounded, a result 
today known as the Bell inequalities. Quantum mechanics predicted that the inequalities 
would be violated, thus giving rise to a way of testing whether the predictions of quantum 
mechanics or the assumptions of Einstein et al. were correct. Experiments have later agreed 
with quantum mechanics |ADR82j . although there are still critics arguing that there are 
loopholes in the experimental assumptions |THf)2| [Sann4j . 

Entanglement has since been regarded as a real - albeit strange - phenomenon of quan- 
tum mechanics. It was widely regarded as being the same as violation of some Bell inequality. 
However, Werner |Wer89j showed that there exist mixed quantum states that allow for a 
local hidden-variable theory, but nevertheless are entangled. This led to new criteria to 
decide if a state was separable Per96, HHHQSj . 

Parallel to this, the view of entanglement had gradually changed. Prom being regarded 
as a curiosity, it was now seen as an information processing and communication resource 
that could be used for performing tasks that would be impossible without it. Among the 
possibilities introduced were quantum cryptography |Eke91j , dense coding |BW92j , telepor- 
tation of a quantum state |BBC^93| and exponential speed-up of certain computational 
tasks |Sho94l ISho97| inro96[ l(Tro97j . Because entanglement was now viewed as a resource, 
it was natural to try to quantify the entanglement in quantum states. 

After the first papers on the quantification of entanglement |SM95| lBBP"'"96l IBBPS961 
IBDSW96,, the subject has grown into a whole field of research. The aim of this work 
is to describe the entanglement measures that have come out of this. The main focus is 
on entanglement between two parties with finite dimensional Hilbert spaces, as this is the 
situation where most of the theory is known. Extensions to entanglement between more 
than two parties and entanglement with continuous variables are also briefly discussed. 

In the course of writing this, several reviews of entanglement and quantum information 
theory have been of great help. The textbook by Nielsen and Chuang |NC00j gives a nice 
introduction to the prerequisites to understand entanglement theory, but contains little on 
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Chapter 1. Introduction 



entanglement itself. John Preskill's lecture notes |Pre98j explain much of the same, but 
has a whole chapter devoted to entanglement. The tutorial by Brufi |Bru02j gives a short 
introduction to most aspects of entanglement, including some entanglement measures. Both 
the lecture notes by Eisert |Eis03j and the review article by Keyl |Key02 include a chapter 
on entanglement measures, but the latter is rather technical in nature. The review article 
by M. Horodecki about entanglement measures [HorOlj summarises the most of the field in 
an elegant way, but it requires much prior knowledge from the reader, and treats almost 
exclusively two-party entanglement. 

The rest of this work is organised as follows. Chapter [2 gives an introduction to the 
concepts needed in the following chapters. A reader with a basic knowledge in quantum 
information theory may skip large portions of this chapter. Chapter describes two-party 
entanglement qualitatively, preparing the ground for chapter |3] which describes the ways 
to quantify two-party entanglement. Finally, in chapter [21 it is briefly indicated how the 
notions from the previous chapters may be extended to systems of continuous variables and 
to multi-party systems. 

Unless stated otherwise all logarithms are taken base-2. 



Chapter 2 



Preliminaries 



Before we can understand how to quantify entanglement, we need to have some concepts 
and procedures clear. In this chapter we will have a first look at qubits and entangled states, 
get used to density operators and see how quantum states can be transformed into others. 

2.1 The qubit 

The simplest nontrivial quantum mechanical system is a two-state system which can be 
described by a vector in two-dimensional complex Hilbert space. Such a system is called a 
qubit. The favourite qubit of many physicists is the spin of a spin-1/2 particle, for instance 
an electron. We will use spins of such particles to illustrate entanglement. 

If we measure the z-component of the spin, the result is either up or down. This mea- 
surement corresponds to the observable operator which we call Z. After the measurement 
the state of the particle is an eigenstate of the observable. We denote these states by \1z) 
and \iz), for spin parallel to and antiparallel to the z-axis. They are mutually orthogonal, 
and we take them as the unit basis vectors of the spin Hilbert space. In addition to the 
eigenstates, the particle may be in any superposition of the eigenstates, 

\^)=a\U)+P\iz) (2.1) 

such that |ap + = 1. When the spin along the z-axis is measured, we obtain +1 with 
a probability of |ap and —1 with probability (where the spin is in units of h/2). After 
the measurement the state is collapses into the corresponding eigenstate. 

As the overall phase of a quantum state does not have any physical significance, only 
the phase difference between a and (3 is important, and for convenience we often choose a 
to be real. In order to also include the property |ap + = 1, we can write the spin state 
(|2.1|) of the particle as 

|^)=cos^|T,)+ei^sin^|U). (2.2) 

This gives a nice and intuitive picture of the state space of the spin-1/2 particle, and thereby 
any qubit. All possible states are covered when the parameters are restricted to < < 27r 
and < 6 < IT. Thus, any state corresponds to a point on the unit sphere, with the 
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azimuthal angle given by (f) and the polar angle given by 9. Such a sphere is called a Block 
sphere, and is a very useful tool for visualising states of a qubit. A Bloch sphere with some 
states indicated is shown in figure 12.11 Note that for 6 = tt, the phase of ()2.2|) changes 
as (j) changes, but the state is still the same. This is reflected on the Bloch sphere, as the 
azimuthal angle is irrelevant at the poles. 




U.> = 75(IT.)-i|i.» 



7;(IT.>+i|W) 



Figure 2.1: The Bloch sphere of the Hilbert space spanned by jf^) and |J,2 



We may choose another set of orthonormal basis vectors. For instance 

ITx) = ^(IU + IU)) IU> = ^(|T.>-|U)). 



(2.3) 



These are the eigenvectors of the X-operator, and they are the resulting states if we measure 
the spin along the x-axis. Likewise, we may measure the spin along any axis, and find a 
corresponding set of eigenvectors. If we denote the basis vectors 



(2.4) 



then the operators for the spin along the x, y and z axes are nothing more than the Pauli 
matrices 



X 



(2.5) 



For the spin along any other unit vector, say ax + by + cz where + b'^ + (? = 1, the spin 
operator is aX + hY + cZ. All operators of this form have eigenvalues ±1 corresponding to 
spin parallel and antiparallel to the vector. 



2.2. Entangled states 
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As already mentioned, other systems are equally good qubits as a spin-1/2 particle. 
Among the more popular ones are the polarisation of a photon and the energy levels of a 
two- level system. They can be represented by a Bloch sphere just as well as the spin of an 
electron can. The advantage of thinking about a spin-1/2 particle is that the vectors on the 
Bloch sphere corresponds to actual spin directions in real space. Thus, a qubit in a state 
close to the y-axis on the Bloch sphere, has, when the qubit is a spin-1/2 particle, a spin 
pointing in a direction close to the y-axis. 

In general the particular implementation of the qubit is not important for our purpose, 
so we simply denote the basis vectors |0) and |1), which in the case of a spin-1/2 particle 
are typically \1z) and \iz), respectively. The superpositions of these states corresponding 
to the spin along the x-axis, are often denoted |-|-) and |— ), according to the sign in the 
superposition. 

|+) = i=(|0) + |l)) |-) = -L(|o)-|l)) (2.6) 

Three qubit operators are of special importance, namely the unitary Pauli operators of 
(|2.5|) . In the context of general qubits, they can be written as 

|1)(0| + |0)(1|, (2.7a) 
y^i|l)(0|-i|0)(l|, (2.7b) 
Z^|0)(0|-|1)(1|. (2.7c) 

X is also called the bit flip operator as it turns a |0) into a |1) and vice versa. Z is called 
the phase flip operator as it switches between the states a |0) -|- /? |1) and a\0) — (3\1). The 
Pauli operators are defined with respect to a particular basis. For instance, in the {|+) , |— )} 
basis, Z takes the role of the bit flip operator. 

The name qubit comes from quantum bit and it plays the same role in quantum infor- 
mation theory as the bit does in classical information theory. The bit is a system which 
can be in one of two states, usually denoted and 1. Thus, each bit is capable of storing 
enough information to answer exactly one yes/no question. As we have just seen, the qubit 
is a much richer structure. In addition to being capable of being in the states |0) and |1), 
just as a bit, it can be in a continuum of superpositions between those states, given by two 
parameters. Still, only a single bit of information can be read from a single qubit. When 
the spin of a spin-1/2 particle is measured along an axis the only possible results are parallel 
and antiparallel to the chosen axis. After the measurement the state is changed, and from 
the result of the measurement we know what it is. No more information about the original 
state can be retrieved. 



2.2 Entangled states 

From the simplest quantum systems whatsoever, we now turn to look at the simplest com- 
posite quantum systems. A composite system is a system which consists of two or more 
parts, and the simplest one is a system consisting of two qubits. We call the two systems 
A and B. In communication and entanglement theory the convention is to think of the 
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subsystems as being in the possession of two well separated observers, called Alice and Bob. 
Any state of each of the two systems can be written as 

\i^)A = « |0)a + |1)a \<t>)B = 7 |0)b + 'J (2-8) 

with |ap + = 1 and |7p + |(5p = 1 The composite state of the two systems is then 
simply the tensor product (or direct product) of the two states. 

I^prod) = (2-9) 

Such a state is called a product state, but product states are not the only physically realisable 
states. If we let the two systems interact with each other, any superposition of product states 
is realisable. Hence, a general composite state can be written as 

\'^)=Y.'^^J\i^^)A®\^J)B (2-10) 

where ^ \otij\'^ = 1 and the sets and are orthonormal bases for the two sub- 

systems. Any composite state that is not a product state is called an entangled state. 

A composite quantum state consisting of two parts only, is called a bipartite state, as 
opposed to multipartite states which consist of more than two parts. For bipartite qubit 
states, four entangled states play a major role, namely the singlet state 

I*-) ^ i=(|01) - |10)) (2.11a) 



V2' 



and the three triplet states 



|^&+) = i=(|01) + |10)) (2.11b) 
1$-) ^ i=(|00) - 111)) (2.11c) 
|$+) = -^(|00) + |ll)) (2.11d) 

where we have used \ij) as a shorthand notation for |i) (g) |j). They are called Bell states 
or EPR pairs. Together they form an orthonormal basis for the state space of two qubits, 
called the Bell basis. The Bell states are maximally entangled and one can be converted 
into another by applying a unitary transform locally on any one of the subsystems. 

Note that if we measure the state of one qubit in a Bell state (that is, measure the Z 
operator which has eigenvalues ±1), we immediately know the state of the other particle. In 
the singlet Bell state, a measurement of qubit A will yield one of the eigenstates |0) and |1), 
each with a probability of 1/2. These results leave qubit B in state |1) or |0), respectively. 

For a single qubit we could always change to another basis where the outcome of a Z 
measurement would be given. For a spin-1/2 particle this means that the spin is always 
pointing in some direction, even though the state will show up as a superposition in a 
basis where the state is not one of the basis states. If the particle is entangled with another 
particle, though, the direction of the spin of that particle alone is not well defined. Actually, 
for particles in one of the Bell states, the probability for measuring the spin of the particle 
to "up" (while ignoring the other particle) is 1/2 for any direction. 



2.3. Density operators and mixed states 
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2.3 Density operators and mixed states 

We want to be able to describe the state of a subsystem that is entangled with another 
subsystem to which we do not have access. State vectors cannot be used for this purpose, 
and we need another representation of quantum states, namely density operators. This is 
not only useful for describing "locally" a subsystem which is part of an entangled system for 
which we know the state. The formalism is also necessary for describing quantum mechan- 
ical experiments - where noise is inevitable - and doing quantum statistical mechanics. In 
the laboratory it is impossible to isolate the quantum systems under study completely. The 
systems become entangled with the environment through unwanted, but unavoidable, inter- 
actions. Nevertheless, it is necessary to describe the system without taking the environment 
into account. 

In the density operator formalism, we describe quantum states by operators on the 
system's Hilbert space instead of unit vectors on it. For any quantum state vector 
the corresponding density operator is the projection operator (^|. As linear operators 
may be - and often is - represented by matrices, density operators are also called density 
matrices. 

So far we have only represented the state in a new way, using density operators instead 
of vectors. But what if we do not know exactly in what state our system is? Imagine, 
for example, that we have a machine in our lab which outputs particles in state (think 
of it as some component of a spin). But once in a while events beyond our control and 
knowledge (a spike in the electric grid, the turning on of an electromagnet in a neighbouring 
lab, etc.) make the machine malfunction slightly producing states which are close to, but 
not exactly |^). Let us for simplicity say that the state {ip') is always created when the 
machine malfunctions, and that the probability for that to happen to any given particle is 
p. Given this uncertainty, the state cannot be described by a state vector. Still, in terms of 
density matrices it is described as 

p=il-p)\i;) {i^\+pm 

A state that may be represented by a unit vector is called a pure state. For a pure state 
we have maximal knowledge of the state. Other states are convex com,binations^ of pure 
states, and are called mixed states. To check if a given density operator, p, represents a 
pure state, it is sufficient to check if p = p^, which holds for all pure states and no mixed 
states. 

The density operator does not tell us from what pure states the state was prepared. For 
example, we can have a state |0) or |1), each with probability 1/2. The density operator 
is then ^(|0) (0| + |1) (1|) = ^I. The same density operator can be made by mixing the 
states 1+) = :^(|0) + |1)) and |— ) = :^(|0) — |1)) in equal amounts, or the three states |0), 

1+) and ^ |0) — |1) in amounts Pi = ^ — P2 = P3 = 1 ~ respectively. 

We call these different collections of pure states with corresponding probabilities, {pi, \ipi)}, 
ensembles. Sometimes it is also called a realisation and we say that an ensemble realises a 

^ A convex combination of elements xi,X2, ■ ■ ■ ,X]\r is an element which can be written ^^PiXi for Pi > 
and ^^Pi = 1. A convex combination of two points is a point on the straight line connecting them. 



8 



Chapter 2. Preliminaries 



mixed state. Hence, when we describe a state from an ensemble by its density operator, we 
discard the information about which ensemble the mixed state was made from. The density 
operator still describes the mixed state as well as can be done, as states from different 
ensembles having the same density operator are experimentally indistinguishable. 

This indistinguishability makes the density operator a perhaps more intuitive represen- 
tation of a state than the state vector. For example, a state vector has an arbitrary overall 
phase factor, so two vectors differing by a phase factor represent the same state. Still, 
when the vector is put together with its dual vector to form a pure state density opera- 
tor, the phase factor vanishes because of the complex conjugation. Hence, all state vectors 
representing the same state are represented by the same density operator. 

To complete our discussion about density operators, we should take note of some of 
their properties. For a more thorough discussion, see e.g. |NC00j . Density operators are (i) 
positive'^, meaning that for any vector |t;), {v\ p\v) > 0, or equivalently that it is Hermitian 
with nonnegative eigenvalues. Also, (ii) the trace of a density operator is unity, tr(p) = 1. 
It turns out that any operator satisfying these two criteria can be realised by a pure state 
ensemble. Hence, we may take these criteria as the defining properties of a density operator. 

We will now see how the density operators are used to describe individual parts of an 
entangled system. In fact, we will see that even when the composite system is in a pure 
state, if the subsystems are entangled, the individual subsystems are in mixed states. 

Consider a bipartite composite system of two qubits, shared between Alice and Bob. 
The qubits are entangled, and the system is in the state 



Now, imagine that Alice loses contact with Bob, and wants to describe her system without 
any reference to the qubit in Bob's lab. Whatever Bob happens to be doing to his qubit 
should not have observable consequences in Alice's lab. If it had, it could have been used 
to perform superluminal communication in some reference frame which would wreak havoc 
in physics by allowing consequences to happen before the cause^. One of the things Bob 
might choose to do with his system, is to measure Z in the {|0) , |1)} basis. If he does, he 
will get |0) or |1) with probability |ap or respectively. This leaves Alice's qubit in the 
same state as the one Bob measured, but as she wouldn't know the result, the state of her 
qubit would be a mixture of |0) and namely 



This mixture correctly predicts the probabilities for all measurements that can be performed 
locally by Alice, regardless of what Bob really does to his system. 

Formally, ignoring some degrees of freedom in a system (like the ones corresponding to a 
particle out of reach) is done by tracing out the relevant degrees of freedom from the density 

^Positive operators here refer to positive semidefinite operators. Positive definite operators, on the other 
hand, satisfy {v\ p > for any \v) and have only positive eigenvalues. 

^This is of course in itself no compelling reason. Even John Bell was afraid that entanglement might do 
away with relativity IWhi98l . The real reason is that quantum mechanics tells us that the probabilities for 
outcomes of local measurements on one system do not change when an arbitrary operation is performed on 
the other system. 






(2.13) 



2.3. Density operators and mixed states 
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operator. This is also called taking the partial trace with respect to the degrees of freedom 
that are to be traced out. Any bipartite state can be written as 

/O = ^ Cijki \ai) {aj\ (g) \bk) {bi\ 
ijkl 

and the partial trace with respect to system i? of a bipartite system is defined as 

trsip) = ^QjfcZ \ai) ® tr ( \bk) {bi\ ) 

ijkl 

= '^Cijki {bi\bk) \ai) {aj\ 

ijkl 

= ^Q,|ai)(a,| (2.14) 

with Cij = Ylki ^ijki {bi\bk)- The density operator obtained by tracing out one part of the 
system is called the reduced density operator. 

We can perform the partial trace on our example state (|2.12j) . The density matrix is 

p = lap |00) (00| + a/3* |00) (11| + a*P\n) (00| + |/3p |11) (11| , (2.15) 

so the reduced density matrix becomes 

PA = tTsip) 

= lap (0|0) |0) (0| + a(3* (0|1) |0) (1| + a* (3 (1|0) |1) (0| + (1|1) |1) (1| 

= |a|2|0)(0| + |/?|2|l)(l|, (2.16) 

which is just what we expected. 

In section [2.21 we saw what we meant by an entangled state as long as the states were 
pure. But this leaves the question of when a mixed state is entangled. Consider for instance 
the mixture of two Bell states 

p = ^\^+){^+\ + ^\^-){^-\. (2.17) 

Written out in an orthonormal product basis, this is 

p = i(|00) (00| + |00) (111 + 111) (00| + 111) (111) 
+ 1(|00) (001-100) (111 -111) (001 + 111) (111) 
= ^|00)(00| + ^|11)(11| (2.18) 

which is a mixture of the product states |00) and |11). So the mixed state can be realised 
by both an ensemble of maximally entangled states and an ensemble of product states. We 
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say that a mixed state is separable if and only if it can be realised as a mixture (convex 
combination) of locally prepared states^. That is, it can be written as^ 

Y^PiP^^^pf^ (2.19) 

i 

where {pi} forms a probability distribution; Pi > and YliPi ~ ^- This also means 
f (HHH96) l that the state can be written as a mixture of pure product states, 

(2.20) 

kl 

Hence, a mixture of entangled states need not be entangled, but a mixture of separable states 
is always separable. Mixture is a process which destroys entanglement. This is because by 
discarding information about which of a number of entangled states the system is in, it can 
no longer be distinguished from a mixture of separable states. 



2.4 The Schmidt decomposition 

A useful tool when working with bipartite, pure state entanglement is the Schmidt decom- 
position. It is a decomposition into the biorthogonal basis which gives the smallest possible 
number of terms for a product basis. 

Given a bipartite pure state \ip), we may write it in terms of some product basis, or- 
thonormal in both subsystems 

i=i j=i 

The coefficients {cy} may be seen as elements in a x dB matrix C. Using the singular 
value decomposition this may be written as 

C = UDV (2.22) 

where D is a dA x ds matrix which is zero except for the diagonal elements which are real 
and positive, U is a dA x dA unitary matrix and V a dB x ds unitary matrix. This means 
that Cij may be written as ^™'^('^'4.rfs) Ui^dk^Vkji so the state vector may be written 

dA dB min(dA.rfs) 

IV') = ^ ^ ^ UikdkkVkj ki) l^i) 
i=l j=l k=l 

min{dA,dB) / dA \ / dg 

= ^ 4fc I ^Uik \ai) I "^Vkj \bj) 

k=l \i=l J \ 1=1 



^This definition is due to Werner |Wer89j who called them "classically correlated" states. 

''In fact it was considered sufficient that that such a sum approximated the state arbitrarily well. It was 
shown in |HHH96) . though, that in the finite dimensional case any separable state could be expressed in cf 
terms, where d is the dimension of the Hilbert space of the composite system. 



2.5. Generalised measurements 
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We now take as our new basis vectors 




) 




and define = dkk to obtain the final form 



min{dA.,dB) 



(2.23) 



fe=i 



which is the Schmidt decomposition. What we have achieved is to write any state as a 
linear combination of maximum min(dyi, ds) product vectors. Also, not only are the product 
vectors in the basis orthogonal, the basis vectors from each subsystem are only used once. 
These vectors are still orthogonal as we have only performed a unitary transformation on 
the set of basis vectors, making the set of new basis vectors biorthogonal. The coefficients 
{Afc} are called the Schmidt coefficients, although in the literature this name is sometimes 
given to {A^} 

For many purposes the mere existence of the Schmidt decomposition is enough to make 
it usable, whereas it sometimes needs to be calculated explicitly. An example of the former 
is to show that the eigenvalues of the two reduced density matrices of a pure, bipartite state 
are equal. Actually, it can be seen directly from (|2.23j) that by forming the density matrix 
and tracing out any one of the two subsystems, the reduced matrix is already on diagonal 
form with {A^} as eigenvalues. It is often easier to calculate the Schmidt coefficients this 
way than doing the singular value decomposition explicitly. 

2.5 Generalised measurements 

To facilitate the quantum operations we will generalise the concept of measurements from the 
way it is usually treated in courses of introductory quantum mechanics. In that formalism, 
which is called projective measurements, the objects representing the measurements are 
Hermitian operators. The result of the measurement is an eigenvalue of the observable and 
after the measurement, the new state of the system is in the eigenstate corresponding to the 
eigenvalue (or projected into the eigenspace if the eigenvalue is degenerate). The probability 
for an outcome of the measurement is the absolute squared expansion coefficient when the 
state is expanded in the eigenstates of the observable. 

In the general measurement scheme, a measurement is represented by a set of linear 
operators {Mi} called measurement operators satisfying the completeness relation 




(2.24) 



m 



The probability for the outcome to be m is 
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and after the measurement with result m, the new state of the system is 



(2.26) 



We may, at least in theory, implement any measurement that can be described by a set 
of measurement operators. Thus, we may construct measurements which convert an initial 
state to a desired state with a certain probability, discarding the system if that particular 
outcome does not occur. As we will see, this is one of the ideas behind distillation of 
entanglement. 

The projective measurement formalism is a special case of the generalised measurement 
scheme. The observable, O, is a Hermitian operator and can be decomposed according to 
the spectral theorem. 

O = '^Xm\Vm) {Vm\ (2.27) 
m 

where {|wm)} are the eigenvectors of the observable, and {Am} the corresponding eigenvalues. 
By taking as measurement operators Mm = |wm) {vml, the completeness relation ()2.24() is 
satisfied. The probability for a given outcome, m, of a measurement is 

{tjj] Ml^Mm IV") = (^bm) {Vm\Vm) (t^mlV') = I ('WmlV') 1^ 

which is exactly the absolute square of the expansion coefficient of for that eigenstate. 
The state after the measurement is 

Mm\'4') _ {Vm\'<lj) \Vm) _ i0 , v 

I — e \Vml 



t A/r U/A I (^'^ 



MlMm 



with an arbitrary phase 9. Like \vm)) this is an eigenstate of O. 

The generalised measurement scheme can easily be further generalised to density oper- 
ators. For pure states, the generalisation is straightforward. The density operator corre- 
sponding to the pure state vector is p = \^) {^\. The probability for a given outcome is 
of course still the same as in 1)2. 25() . but we write it in terms of the density operator: 

p{m) = tv{MmpMl) (2.28) 

When the state changes after a measurement with result m, the density operator corre- 
sponding to the state vector (|2.26l) is obviously 

Mm Mi ^ MmpMl ^2 

(V'l MlMm IV') tliMmpMjn) 

The probability (|2.28j) and new state (|2.29|) is also valid if the density operator represents 
a mixed state. To see this, we write the density operator in terms of one of its pure state 
realisations: 
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The joint probability for the system to be in the pure state \tpi) and measuring a result m 
is pip{m\i). Summing over all i we get the total probability for getting m as the result of 
the measurement. 



The general measurement scheme covers all operations that can be performed on a quan- 
tum system. A unitary transformation is the special case of a one-outcome measurement as 
the measurement operators then will have to be unitary to satisfy the completeness relation 
1)2. 24() . Performing a different unitary transformation Um conditional on the outcome of 
a measurement described by the operators {Mm}, is described by the new measurement 
operators {UmMm}- It is easy to see that these operators satisfy the completeness rela- 
tion, that the probabilities of the outcomes are the same and that the post measurement 
state is what we expect, namely Um{MjnpMm)Um- Similarly, consecutive measurements 
can be put together to form a single set of measurement operators, even when the mea- 
surements to be done depend on the outcome of the previous ones. When doing this, all 
possible combinations of measurement outcomes are considered as single outcomes of the 
total measurement. 

In applications it depends on the situation whether we use the full generalised mea- 
surement scheme or only consider projective measurements. Projective measurements are 
sometimes simpler to handle and combined with unitary transformations they describe ev- 
erything that can be described by generalised measurements. Therefore, both formalisms 
will be used in the following. 

2.6 Transformation of states and LOCC 

We will see later that the set of operations that we allow to be performed on a quantum 
state is essential to the most basic entanglement measures. All operations can be described 
as generalised measurements. But we will find it useful to describe them as measurements 
combined with unitary transformations where the next step in a sequence of operations may 
be decided by a previous measurement outcome. 

2.6.1 Local Operations and Classical Communication 

In the context of entanglement a subclass of quantum operations is particularly important, 
namely local quantum operations (LO). When we have a bipartite or multipartite state, the 




ti{MmpMl) 



Similarly, the new state after the measurement of the mixed state becomes 
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subsystems may be separated, or at least well isolated from each other. We may therefore 
only perform operations that act locally on each subsystem. That is, the operators will 
be of the form Oa ® Ob ® ■ ■ ■ ■ In addition, the parties sharing a multipartite state are 
usually allowed to communicate classically, and through this classical communication (CC) 
they may perform different operations depending on the measurement results of the other 
parties. This subclass of operations is called Local Operations and Classical Communication, 
usually shortened to LOCC. In the literature it is sometimes named otherwise, e.g. QLCC 
or LQCC. LOCC is important in entanglement theory because entanglement is exactly those 
correlations that cannot be created using LOCC. 

2.6.2 Transformation of a Bell state 

Suppose Alice and Bob share a Bell state |$^) = |00) + |11) and want to transform 
it into the state cos 6 |00) + sin 6 |11). One way of doing this is as follows. Alice performs a 
measurement described by the measurement operators 




(2.31) 



where the local state vectors are denoted as usual: 

10) = (;) ID = (?) . (2.32) 

The operators clearly satisfies the completeness relation (|2.24j) . as 
/cos2 e \ /sin2 e \ _ A 0\ 

V sin2 e) ^ \ cos2 e) ~ \0 l) ' 

After the measurement we may have one of two states, depending on the outcome, 

iV^i) = cos6'|00) + sin6l|ll) IV'2) = sin6l |00) + cos6l |11) 

where is the desired output state. If the result of the measurement is 2, leaving the 
system in the state \ip2), Alice applies the X operator (also called bit flip operator) to her 
system converting |0) into |1) and vice versa. Then she tells Bob to do the same through 
the classical communication channel. The resulting state is and their common goal is 
achieved. 

The transformation we just did, transforms a maximally entangled state into a less or 
equally entangled state. For 6 = and = 7r/2 the output state is a product state, which is 
not entangled at all, whereas for = 7r/4 and 9 = Sir /A the output state is also a Bell state. 

2.6.3 Concentrating partial pure state entanglement 

We may also try to do the transformation the other way around. We begin with the state 
where we for simplicity suppose that < < 7r/4, and want to end up with the Bell 
state l^*"*"). The lower boundary of this interval is a product state for 9 = which cannot 
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possibly be turned into an entangled state by LOCC. The upper boundary = 7r/4 is the 
desired output state and will need no transformation. 

The problem is now to construct local measurement operators which give \^~^) as a post 
measurement state. The operator 

Ml = const • J/] 

will give the desired output, and we want the constant to be as large as possible to maximise 
the probability for this output. Still, for the completeness relation (|2.24() to be satisfied none 
of the diagonal elements can be greater than 1. This is because for any Mm, MmMm has 
nonnegative diagonal elements, and each diagonal element needs to sum up to 1 to get the 
identity operator in (|2.24l) . In the valid range of 0, cos 9 > sin 9, so the maximal constant is 
sin 9. This gives the optimal measurement operator for the successful outcome 

^^=(^T ?) =tan0|O)(O| + |l)(l|. 

The remaining set of matrices representing the measurement operators need to have a van- 
ishing element to the lower right, and the remaining elements need to be such that the 
completeness relation is satisfied. The easiest way to do this is to simply add one more 
operator, namely 

The new states after the measurement are proportional to 

Ml iV'i) = sin6l|00) + sin6l|ll) 
M2 iV'i) = cos6l\/l -tan2 6'|00) . 
The probability for the two outcomes are 

= {ipil mJMi iV'i) = sin2((00|00) + (11|11)) = 2sin2 6l 
p{2) = {ipi\ m|M2 iV'i) = cos^ 61(1 - tan^ 9) = cos^ 9 - sin^ 9 = 1 - 2sin2 9 

and the normalised new states are 

|V'out,i) = ^(|00) + 111)) = !<!>+) 

|^out,2) = |00) . 

We see that there is a certain probability for converting a partially entangled state into 
a maximally entangled state, but the less entangled the initial state is, the smaller is the 
probability. This is consistent with the fact that LOCC cannot increase the expected entan- 
glement. Still, we may perform a measurement which may give us a more entangled state at 
the risk of losing the entanglement we had. If we have a large number of bipartite two-qubit 
systems, we may apply the above procedure and discard the systems where the output is a 
product state, thereby effectively concentrating entanglement. 
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2.6.4 Distilling entanglement from mixed states 

Also mixed state entanglement may be concentrated, or distilled. This will be important for 
practical applications, as the channels used to create entanglement in general will be noisy. 



A collection of distillation protocols was introduced in BBP^96 and |BDSW96| . We will 
take a closer look at a slightly modified version of a distillation protocol introduced in the 
latter, which for certain mixed states can yield pure Bell states as output. 

Consider a mixed state, as always shared between Alice and Bob, which can be rep- 
resented as a collection of l^*"*") states contaminated with a portion of |01) states. The 
corresponding density operator is then 

p={l-p) |01) (01| («>+| 

= (1 - p) |oi) (Oil + 1 |oo) (oo| + 1 |oo) (111 + 1 111) (oo| + 1 111) (111 . 

Alice and Bob now pick two pairs from the ensemble. Each of them will perform a joint 
operation on the two subsystems on their side. Then they measure one of the pairs, thereby 
destroying any entanglement in that pair, to decide if the other pair is in the desired output 
state. 

The composite system of the two pairs is described by the density operator 
p®^ = 

(1 -p)'' 10101) (01011 + 10100) (01001 + £^ 10100) (01111 + £^ join) (01001 + £^ 10111) (01111 

+ ^^^2^ 10001) (0001| + ^ loooo) (ooool + ^ loooo) (0011| + e1 lOOll) (ooool + ^ lOOll) (0011| 
+ ^^^2-^ lOOOl) (1101| + ^ loooo) (llOOl + ^ loooo) (1111| + ^ lOOll) (llOOl + ^ lOOll) (1111 
+ ^ii^^ |1101) (0001| + ^ |1100) (OOOOl + ^ |1100) (0011| + e1 |llll) (ooool + el |iiii) (0011 
+ £(!_£) 11101) (11011 + si 11100) (11001 + si 11100) (11111 + si 11111) (11001 + e1 |1111) (11111 

where the first and third qubit are in Alice's possession and the second and fourth are in 
Bob's possession. The two first qubits are the first pair, and the two last are the second 
pair. It may be written symbolically as \ A1B1A2B2) ■ Both Alice and Bob applies a CNOT 
gate to the two qubits in their possession, using the qubit from the first pair as a control 
bit and the second one as the target. A CNOT gate is a unitary transformation which does 
nothing when the control bit is zero and flips the target bit if the control bit is set to 1. The 
operator is written 

f/cNOT = |00) (00| + |01) (01| + |10) (11| + 111) (10| 
= |0) (0| 01+ |1) (1| 0X 

where X is the usual bit flip operator, or Pauli X-operator. But as Alice's operator works 
on the first and third qubit while Bob's operator works on the second and fourth, in our 
basis the operators are 

C/cnot(1,3) = |0)(0|® I + |1)(1|® I (g)X0l 

C/cnot(2,4) = I ®|0)(0|®I®I + I «)|1)(1|®I®X 
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Operationally it only amounts to flipping the third qubit of p®^ in all terms where the first 
qubit is 1 and flipping the fourth qubit where the second qubit is 1. This leaves the new 
state 



Pout = Z7cNOT(l,3)[7cNOT(2,4)p®2c7t^Q^(2,4)[7tj,oT,(l,3) 



(l-p)^|0100) (0100| + 



p(i-p) 



|0101) (0101| + 



p(l-p) 



|0101) (0110| + 



p(l-p) 



|0110> {0101| + 



p(l-p) 



|0110> (0110| 



+ 



+ 



+ 



p(l-p) 

2 

p(l-p) 



p(l-p) 



|0001) (0001| + ^ |0000) (ooool + ^ lOOOO) (0011| + ^ |0011) (0000| + 



^ |0011)(0011| 



lOOOl) (lllOl + ^ 10000) (1111| + ^ loooo) (iiool + 



10011) (Hill 



+ 2^ 10011) (11001 



|1110) (0001| + ^ |1111) (0000| + 



^11111) (00111 



+ ^ |1100) (0000| + ^ |1100) (0011| 



+ 



p(l-p) 



|1110) (1110| + 



^11111) (11111 



+ ^ |1111) (1100| + \ |1100) (1111| + ^ |1100) (1100| . 



In this state Alice and Bob both measure the state of their qubits belonging to the target 
pair. This can be done formally with the set of measurement operators 



Moo = I®I®|0) (0|®|0) (0| 
Moi =I®I®|0) (O|0|l)(l| 
Mio = l0l0|l)(l|0|O) (0| 
Mil = 1(8)1® |1) (1| |1) (1| 



(2.33a) 
(2.33b) 
(2.33c) 
(2.33d) 



and is equivalent to measuring the observable I (X" I (g) ^3 ^4 as the measurement operators 
are constructed from eigenvectors of this observable. The interesting outcome is when both 
qubits are measured to 1. The probability for that is 



p{ll) = tr (MiipoutM^i) 



(2.34) 



where only the boxed terms of pout survive being operated on by the Mn and of them 
only the two diagonal ones give a nonzero contribution when the trace is taken, giving a 

2 

probability of ^. The new state after the measurement is 



= |$+) ($+1 (g) |11) (11| . 



(2.35) 



So by using two entangled pairs in a mixed state Alice and Bob have obtained a pure 

2 

state where one of the pairs is in a Bell state. The probability for this to happen is so 

2 

on average we get a yield of ^ Bell states per input pair. 

Some additional Bell states may be squeezed out of the mixed states after using the 
above protocol. The measurement result 00, with a probability — 2p + 1, leaves the pair 
of control qubits in the state 

2(1 -p)2 



Poo 



3p2 -4p + 2 



|01) (01| + 



3p2 _ 4p + 2 



|$+) 



(2.36) 



18 



Chapter 2. Preliminaries 



which is on the same form as the original input state with a new mixing ratio 



P 



Sp2 - 4p + 2 



(2.37) 



These pairs may be run through the same procedure once more with a yield of increasing 

the total yield with an amount of 2^(00) " where the factor ^ comes from the fact that 
we use two pairs and only get one out. This procedure may be iterated, and the total yield 
converges quickly so only few iterations are necessary in practice. Figure 1221 shows the total 
yield calculated numerically as a function of p. 



0.35 




Figure 2.2: The yield of the distillation protocol for 1, 2 and 3 iterations. The dotted line is the 
asymptotic limit. 



In general, distillation protocols do not yield pure state pairs after a finite number of 
iterations. Likewise, in the above protocol this only happens when the state can be written 
as a mixture of a Bell state and another special state. A general feature, though, is the 
application of the bilateral CNOT unitary transformation on two pairs. 

The purity of the output states is measured by the fidelity of the state p, defined as^ 



(2.38) 



where |<1>+) is the desired output state. Applying a distillation protocol to an ensemble of 
states, leaves a subensemble where the fidelity is higher, along with a subensemble where the 
fidelity is lower, often zero. In the limit of reapplying the procedure to the higher-fidelity 
states an infinite number of times, the fidelity approaches unity. This comes at a cost. 



Sometimes the fidelity is defined as the square root of this. 



2.6. Transformation of states and LOCC 



19 



however, as the yield decreases as we demand higher fidehty. For the simplest protocols 
the yield approaches zero as the fidelity approaches unity. But by combining protocols and 
using the already distilled states as a resource, the limit can be made finite. 
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Chapter 3 



Characterisation of bipartite 
entanglement 



Before moving on to quantify the entanglement in a bipartite state, we will try to characterise 
it more qualitatively. There are a number of questions which does not have an immediate 
answer. Given a state, how can we tell wheter it is entangled? Is there more than one type 
of entanglement? Can all entangled states be distilled to Bell states? 

3.1 Pure states 

When we are given a bipartite state the first question of interest is whether it is entangled or 
not. This is the question of separability. A separable state is a state which is not entangled, 
so each of the subsystems can be given a description on its own. When it comes to pure 
states, a state is separable if and only if it is a product state. As indicated in 12. 2[ a pure 
state is a product state if it can be written as 



for any pure states \tp)j^ and belonging to Alice and Bob, respectively. If the state is 
not a product state, it is an entangled state. 

In practice, though, one needs to find the right basis vectors to be able to write a product 
state as in 1)3. A simpler criterion is desired. The Schmidt decomposition (section 12. 4|) 
is helpful in this respect. The number of terms in the Schmidt decomposition is called the 
Schmidt number, sometimes also Schmidt rank. That is, the number n in 



The Schmidt number is at most mm(dA, ds), where dA {ds) is the dimension of subsystem 
A {B). But when some of the Schmidt coefficients are zero, those terms may be skipped 
giving a lower Schmidt number. Remember that the Schmidt coefficients are simply the 
square root of the eigenvalues of any of the reduced density matrices. Hence, the Schmidt 
number is simply the number of nonzero eigenvalues pA or ps- 





n 





k=l 
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What the Schmidt number tells us physically is basically how many degrees of freedom 
that are entangled. And as a product state can be written as 1)3. If) it is obvious that 
its Schmidt number is 1. For pure states this is a necessary and sufficient criterion for 
separability. 

It turns out that entanglement from any pure entangled state can be concentrated into 
Bell states. This may be done by using the procedure similar to the one described in 
section 1!^. 6. 31 or a more advanced procedure that operates on more partially entangled pairs 
at the same time to give a higher yield (see e.g. pBPS96 ). Bell states may again be 
used to produce other entangled states by means of LOCC. Because of this, and the easy 
identification through the Schmidt number, the only classification we make for bipartite 
pure states is in separable and entangled states. 



Some properties of mixed states were discussed in section 12.31 The structure of mixed 
state bipartite entanglement is richer than its pure state equivalent. For instance there are 
entangled states that cannot be distilled to Bell states, and cannot be put to use in any 
of the applications where entanglement is a vital resource. Also, the more difficult task of 
finding out whether a given mixed state is separable or not gives rise to classes of states 
which can be identified by other criteria. 

3.2.1 Separable states 

Recall from section 12.31 that a bipartite state is called separable if it can be realised as a 
mixture of product states. 



with Pi > and YliPi ~ ^- Given a density matrix in some basis the task of finding local 
density matrices to satisfy (|3.3|1 can be enormous. Even for the simplest system which may 
be entangled, a pair of qubits, a density operator would typically look like 



3.2 Mixed states 




(3.3) 




(3.4) 



or in matrix representation in the {|00) , |01) , |10) , |11)} basis, 
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It is not easy to see that this can be written as ()3.3() . Still, it is constructed as a mixture of 
two product states. It can be written as 

(A) ^ (B) , (A) ^ (B) 



where 



It could of course be realised by a variety of other mixtures as well. 

This example shows that we could well use an operational criterion for separability, 
like the Schmidt number for pure states. In general, no easily computable necessary and 
sufficient criterion is known. Still, we have a necessary condition which is easy to compute. 
This is known as the PPT^ -criterion or Peres- Horodecki criterion Per96, IIIIII96 . It 
states that if the state is separable, then the partial transpose of the density operator with 
respect to one subsystem is positive. A transpose of an operator needs to be taken with 
respect to a basis, and the resulting operator depends on the basis. Actually, if is 
the transpose of the operator with respect to one basis, any operator of the form UA^U'^ 
for unitary U satisfying U = (in the same basis) is a transpose of the operator with 
respect to another basis. However, since the different bases are connected by a unitary 
transformation, the eigenvalues are independent of the basis. Therefore, any orthonormal 
basis will do for taking the transposition. More formally, if we choose an orthonormal 
product basis {\viVj)} = {\vi) \vj)} for the state p, the partial transpose p^^ is defined by 
its matrix elements 

Pmfi,nii ~ {'^mVfji] p ^ l^n'L'i/) = Pmu,nfj,- (3-7) 

In the example matrix taking the partial transpose of system B corresponds to trans- 
posing within each of the four 2x2 blocks, whereas switching the upper right block with 
the lower left block corresponds to transposing system A. 

For systems of dimension 2x2 and 2x3 the PPT-criterion is also sufficient for the state 
to be separable |IIHII96j . For larger systems, though, there exist entangled PPT states 
|Hor97j . The PPT-criterion is still so useful that the class of PPT states is an important 
class of states, strictly larger than the class of separable states. Anything that can be proved 
for all PPT states, is automatically true for separable states. 



3.2.2 Free and bound entanglement 

Unlike pure states, not all entangled mixed states can be distilled into Bell states |HHH98j . 
This is an important distinction, because for many applications entanglement must be dis- 
tilled in order to be useful. Entanglement that can be distilled is called free entanglement 
and entanglement that cannot be distilled is called bound entanglement. Bound entangled 
states cannot be written as a convex combination of product states (|3.3() and thus cannot 
be prepared locally by two classically communicating parties, but its entanglement is well 

^Positive Partial Transpose 
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hidden and difficult to detect. It has been shown that for a pair of qubits all entangled 
states can be distilled |HHH97j . 

A connection between free entangled states and PPT states was found in HHH98|. 
It says that no free entangled state is PPT. Hence, all PPT entangled states are bound 
entangled states. This result was actually the first evidence of the existence of bound 
entangled states. 

This raises the question of whether a nonpositive partial transpose (NPPT) state is 
always distillable. The answer to this question is not known, but there are indications that 
NPPT bound entangled states do exist jPSS+Onl iPT^LBnOj . 

It has been proven that bound entanglement is useless for quantum teleportation |HHH99aj . 
and it can be proved for other applications as well. One has found, though, that bound en- 
tanglement is not completely useless after all. For certain tasks it can be activated HHH99b] 
to facilitate distillation of other quantum states. It has even been shown that a small amount 
of PPT bound entanglement as a resource makes it possible to distill entanglement from 
any NPPT state \'\\'()2l)j . This means that if NPPT bound entangled states really ex- 
ist, entanglement can be distilled by putting together two states that are not distillable by 
themselves. 



Chapter 4 

Measures of bipartite entanglement 



Entanglement can be utilised for performing many tasks which are impossible without it, 
or enhance the performance of other tasks. In most of these tasks entanglement is con- 
sumed, so we trade entanglement for something else. It is clear from this that we may treat 
entanglement as a resource, just as energy is a resource needed for certain tasks. This is 
the reason why we want to quantify the entanglement in an entangled state. Given a state 
and a task that consumes entanglement, how much can we achieve? How well can we do? 
It is not obvious that an entangled state that performs a given task better than another 
is will be the better resource for another task. Indeed, this is only true in a very limited 
context. Because of this, we have in general many ways to quantify the nonlocal resources 
(or entanglement) in a quantum state. 

An entanglement measure E{-) is a functional that takes a quantum state of a multi- 
partite system to a nonnegative real number. In this chapter we only consider bipartite 
systems, hence 



where V{TL) is the set of density operators on the Hilbert space Ti = TLa '^'Hb- There are 
in principle two ways to quantify the entanglement in a state. Operational measures are 
based on how well a certain task can be performed, usually compared to Bell states. The 
other way, which gives rise to abstract measures, is to work from a set of natural axioms 
we believe an entanglement measure should satisfy, and look for functionals that satisfy the 
axioms. 

There are also two regimes to consider. In the finite regime we consider the resources 
contained in a single quantum system, whereas in the asymptotic regime we take into account 
an infinite number of systems in the same state and consider the the resources per system. 

4.1 Pure states 

Before treating the general case of mixed states we will consider the simpler case of en- 
tanglement measures of pure states. The treatment is facilitated by the fact that a pure 
state contains no classical correlations between the subsystems, so any correlation present 
must be of quantum nature. As mentioned in section 12. 3( even if an entangled state is 



E : V{n) M" 



+ 
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pure the states of the subsystems - described locally by the reduced density operator - are 
mixed. For pure states, then, the amount of "mixedness" turn out to be a good measure of 
entanglement. 

4.1.1 von Neumann entropy 

Our measure for how mixed a quantum state p is, will be the von Neumann entropy^ defined 
as 



where we take the logarithm base-2 as is the custom in information theory. It is most easily 
calculated from the nonzero eigenvalues Aj of p as 



The von Neumann entropy is often seen as a generalisation of the Shannon entropy 
from classical information theory. Although it is true that Shannon entropy arises as a 
special case in quantum information theory when only orthogonal states are considered, the 
historical infiuence is the other way around, von Neumann introduced his entropy in 1932 
|vN32j and Shannon published his mathematical theory of communication in 1948 |Sha48j . 
In fact, it was von Neumann who suggested that Shannon should name his uncertainty 
function "entropy" as that was already used in statistical physics |TM71[ p. 180]. 

What the von Neumann entropy essentially describes, is the uncertainty in a quantum 
state. It is zero for pure states, and smaller for a mixture of two nonorthogonal states 
than a similar mixture of two orthogonal states. It is not only an analogue to the entropy 
known from thermodynamics, it is the same quantity (up to a constant with the definition 
(1121)). Indeed, it was derived using notions such as the ideal gas law and other laws of 
thermodynamics jvN321 p. 191-202] 

4.1.2 Reduced von Neumann entropy — entropy of entanglement 

For the purpose of quantifying entanglement, the important quantity is the reduced von 
Neumann entropy. This is the von Neumann entropy of the reduced density matrix. We 
saw an example in section 1)2. 3|) that the reduced density operator of a pure entangled 
state represents a mixed state. The von Neumann entropy of either of the subsystems (the 
entropies are equal) is a good measure of entanglement. Because of this, the reduced von 
Neumann entropy is also known under the name entropy of entanglement Ee- 

It is easy to see that the reduced von Neumann entropy is equal for both reduced density 
matrices. The von Neumann entropy only depends on the eigenvalues of the density matrix 
()4.3|) and because of the existence of the Schmidt decomposition, the eigenvalues are equal, 
as was mentioned in section 

We haven't yet pointed out what constitutes a good measure, but the entropy of entan- 
glement has some properties we find natural, (i) It is zero for any product state, (ii) it is 
maximal when the reduced density matrix is completely mixed, e.g. when the subsystems 

^Pages 359-379 in the far worse typeset English translation |vN55| . 
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have no individual properties (for the degrees of freedom considered) and (iii) it is invariant 
under local unitary transformations. The above requirements were introduced in j5M95] . 
but requirements have later been refined and placed on firmer ground for both pure and 
mixed states. 

The entropy of entanglement was first introduced as a measure of entanglement in 



BBP"'"96| . It is an abstract measure in the sense that it satisfies some requirements, and 
it does not have an immediate operational interpretation. We will now take a look at two 
pure state entanglement measures which have a direct operational meaning. 

4.1.3 Entanglement cost and distillable entanglement 

We saw in section l^.6.c{l an example of a concentration of pure state entanglement. We turned 
a partly entangled two-qubit state into a Bell state with a yield depending on how entangled 
the initial state was. There are other protocols which give a higher yield if we have more 
than one pair to operate on (see e.g. BBPSOfT). The distillable entanglement Ed is defined 
as the maximum yield of Bell states that can be obtained, optimised over all possible LOCC 
protocols. The distillable entanglement is also sometimes called entanglement of distillation 

The entanglement cost Ec is the dual to the distillable entanglement. Two separated 
parties cannot prepare an entangled state if they can only communicate classically. But if 
they have some entangled state in some standard form like Bell states, they can convert those 
into the desired entangled state. Thus, the entanglement cost is defined as the minimum 
number of Bell states needed to create a given state by means of LOCC. 

The above definitions were implicitly given for the finite regime. It turns out that in 
this regime E^ and Ec are hard to calculate even though the states are pure. What is clear 
intuitively is that for any state 

Ed < Ec. (4.4) 

If not, one could create entanglement by means of LOCC by converting Bell states to a 
state not satisfying (|4.4|) and converting them back again. 

In the asymptotic limit, the results are known. We define the regularised (or asymptotic) 
versions of distillable entanglement and entanglement cost as 

Eoim = E^im ^ hm (4.5a) 

Ecm) ^ E^im ^ hm ^^i^^. (4.5b) 

n— »oo n 

In this case, it was shown in |BBPS96j that both the distillable entanglement and entangle- 
ment cost is equal to the entropy of entanglement, £'£)(!'(/')) = Ec{\'4')) = Ee{\'4'))- In fact, 
if two states iV'i) and \'il)2) have the same entropy of entanglement, they can be intercon- 
verted with efficiency approaching unity as n ^ oo. If they do not have the same entropy 
of entanglement, can be converted into |'(/'2) with asymptotic yield Ee{\iI>i)) / Ee{\iI^2))- 
It has even been shown that this conversion can be done such that the amount of classical 
communication required per Bell state produced approaches zero |LP99j (it scales as ^/n). 

Despite the fact that the entropy of entanglement did not have an operational interpre- 
tation ah initio., we see that it indeed plays a fundamental role. The fact that the entropy 
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of entanglement can be asymptotically conserved in the conversion of states using LOCC, 
places huge constraints on other candidates for pure state asymptotic entanglement mea- 
sures. A measure would not live up to our expectations if we could increase it using only 
LOCC. If another measure were to give two different values for states with the same en- 
tropy of entanglement, it would be possible increase it by LOCC by converting the state 
with lower value to the one with higher value. The conversion between states with con- 
version rate Ee{\iPi)) / Ee{\iP2)) imposes even stronger constraints, so we would expect any 
asymptotic pure state measure to coincide with the entropy of entanglement. The exact 
conditions for this to happen are discussed in the next section. 

4.1.4 The uniqueness theorem for measures of entanglement 

The fact that the pure state asymptotic entanglement cost and distillable entanglement 
coincide with the entropy of entanglement is a special case of a more general property of 
entanglement measures. The uniqueness theorem for entanglement measures states that any 
pure state measure of entanglement that satisfies certain natural criteria coincides with the 
entropy of entanglement. The criteria are considered to be too strict for the finite regime, 
but all natural asymptotic measures satisfy the conditions. 

The uniqueness theorem has been developed gradually in the literature, in the beginning 
only as an intuitive idea of the same character as given in the previous section |FR97| . 
Later, uniqueness was proved from a set of conditions that measures should satisfy |Vidnni 
IHHHOOj . In |DHR02j the minimal conditions for the theorem to hold were found. The 
conditions come in difi^erent versions, some weaker than others. It turns out that some of 
the stronger conditions are not necessary for the theorem, but are in fact implied by the 
weaker conditions. The strongest versions are listed first. 

(PO) If {tp) is separable (i.e. a product state), then Edip)) = 

(Pla) (Normalisation) For the maximally entangled state in d x d dimensions 



where the two sets and are orthonormal bases for the two subsystems, E{\^'^)) = 

logd, where we as usual take the logarithm base-2. 

(Fib) For the Beh state |$+) = -^(|00) + |11)), E{\<^+)) = 1. The previous condition 
is simply the generalisation of this. 

(P2) For any operation A that can be implemented by means of LOCC and any \ip) such 
that A(|^)) is a pure state, E{A{\ip))) < E{\7p)). 

(P3) (Continuity) Let {|V'n)} and {|0n)} be sequences of pure bipartite states living 
on the sequence of Hilbert spaces {Tin}- For all such sequences such that || {ipn) (V'nl — 
|0n) {'Pn\ 111 — ^ where || • ||i is the trace norm ||^||i = tr(\/]4U), 





(4.6) 




l + log(dim?t:„) 



0. 



(4.7) 
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(P4a) (Weak additivity) For all pure states and n > 1, 



n 



E{m (4.8) 



(P4b) (Asymptotic weak additivity) Given e > and a pure state IV'), there exists an 
integer > such that for all n > A^, 

^M!!)_,<E(i«)<£M!!)+, («) 

n n 
(P4c) (Existence of a regularisation) For all bipartite pure states | tp) , the limit 

E°°{\ij)) = lim ^^^"^^^ ^ (4.10) 

n^oo fi 

exists. 

With these conditions given, the uniqueness theorem for entanglement measures states 
|DHR,n2j that for a functional E on pure states, the following are equivalent: 

(i) E satisfies (Plb), (P2), (P3) and (P4b). 

(ii) E satisfies (PO), (Pla), (P2), (P3) and (P4a). 

(iii) E coincides with the entropy of entanglement E = Ee- 

On the other hand, if E satisfies (PO), (Pla), (P2) and (P3), it automatically satisfies (P4c), 
and E°° = Ee- 

(i) represents the weakest set of conditions. When they are satisfied, the stronger con- 
ditions in (ii) are also satisfied. The connection is that (P4a) and (P4b) turn out to be 
equivalent ((P4a) ^ (P4b) is obvious) and (Plb), (P2) and (P4a) together implies (PO) and 
(Pla). The uniqueness for the asymptotic regime, previously mentioned, is concretised by 
the last sentence. 



4.2 Mixed states 

The measures of pure state bipartite entanglement constitute a very limited subset of the 
measures of generally mixed states. For mixed states there is no unique way to quantify the 
nonlocal quantum resources even in the asymptotic regime. In the pure state asymptotic 
limit, the entropy of entanglement imposed a total order on the set of states. Any state 
could be converted into another state with equal or less entanglement. In the mixed state 
case there are states where neither can be converted into the other, giving only a partial 
order on the set of states. (For a longer discussion on entanglement measures as ordering 
of the states, see |MSV04j .') This makes the structure of mixed state entanglement richer 
than that of pure states. 

4.2.1 Entanglement cost and distillable entanglement for mixed states 

The entanglement cost and distillable entanglement can easily be defined for mixed states in 
the same way as for pure states. The entanglement cost Ec is the minimum number of Bell 
states needed to produce the state by means of LOCC and the distillable entanglement E^ 



30 



Chapter 4. Measures of bipartite entanglement 



is the maximum number of Bell states that can be distilled by an optimal LOCC distillation 
protocol. However, a perfect conversion between mixed states in the finite regime is usually 
not possible, so in the following we will only consider the asymptotic versions Ed and Ec- 
There is a tweak, though, compared to the pure state case. As mentioned in the discussion 
of entanglement distillation fsection I2.6.4|) we do not require the conversion to be perfect. 
The only requirement is that the fidelity approaches unity as n — > oo. 

Note that to produce mixed states from Bell states using perfect LOCC we will need to 
discard some information. This can be done by doing a random local unitary transformation 
and not recording which one it was, or by doing a measurement and forgetting the outcome. 
A third option is to entangle the local system with an ancilla system and discard the ancilla. 
It is not surprising then that the distillable entanglement usually is strictly smaller than the 
entanglement cost. Actually, for any measure E which satisfies certain natural conditions 
for asymptotic measures, the entanglement cost and distillable entanglement provide upper 
and lower bounds |HHHOnp . 



The bound entangled states, discussed in section in.2.2[ is an example where the difference 
between Ec and E^ is highly visible. Ed is obviously zero, but Ec is finite. This illustrates 
that one single measure is not enough to quantify the entanglement resources in a mixed 
quantum state, even in the asymptotic limit. 

4.2.2 Axioms for abstract measures 

The axiomatic approach to quantifying entanglement has been quite successful in finding 
measures that quantify certain aspects of entanglement. During the course of finding the 
right conditions, some conditions have been discarded as unnecessary and others have been 
restricted to a specific regime. The common denominator, though, has always been that 
entanglement, irrespective of how we quantify it, should not increase on average by any 
LOCC operation, i.e. it is monotonic under LOCC. The first attempt was done by Vedral et 
al. in |VPRK97] and was slightly improved by some of the same authors in |VF98| . Vidal 
|Vidnfl | argued that the only absolute requirement was monotonicity under LOCC. This au- 
tomatically implies other properties such as nonchange under local unitary transformations, 
convexity and constancy on separable states. 

The conditions on mixed states discussed here are taken from |DHR.n2| . The previously 
discussed conditions on pure states are specialisations of those. 

(EOa) E{p) = if and only if p is separable. This is a useful property for measures that 
satisfy it, but it is too strict in general. Bound entangled states, for instance, have zero 
distillable entanglement, but are not separable. 

(EOb) E{p) = if /) is separable. 

The normalisation conditions are the same as for pure states. 

(Ela) (Normalisation) For the maximally entangled state uidxd dimensions, I'l'^) {^^\^ 

^In HHHOO the upper bound was Ep, the regularised entanglement of formation. This has been shown 
to be equal to the asymptotic entanglement cost |HHTOl |. 



Ed ^ E ^ Ec 
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where 

d 

\^d) = T.^\')A^\i)B (4-12) 

and the two sets and are orthonormal bases, E'(|$^) ($^|) = logd. 

(Elb) For a BeU state |$+) = ^(|00) + |11)), E{\^+) ($+|) = 1. 

(E2a) (LOCC monotonicity) For any LOCC operation A, E{A{p)) < E{p). 

(E2b) When A is a strictly local operation which is either unitary or adds extra dimen- 
sions, then E{A{p)) = E{p). 

(E2c) When A is a strictly local unitary operation, then £'(A(p)) = E{p). 

(E3a) (Continuity) Let {pn} and {(Jn} be sequences of bipartite states living on the 
sequence of Hilbert spaces {Hn}- F or al l such sequences such that ||pn — Cn||i — where 
II • 111 is the trace norm |m|i = A'^A), 

^f"") - ^f""' .0. (4.13) 



1 + log(dim'H„ 

(E3b) The previous condition (E3a) is weakened by only requiring it to be satisfied for 

approximations to pure sates. That is, (E3a) only needs to hold when the p„ arc pure. 

Many measures satisfy the weak additivity condition, but it is not considered necessary 
for a good measure. 

(E4a) (Weak additivity) For all states p and n > 1, 

^-^^Ei,). (4.14) 
n 

(E4b) (Asymptotic weak additivity) Given e > and a state p, there exists an integer 
N >Q such that for all n>N, 

E(£::l.,<E(,)<^ + e. (4.15) 

n n 

(E5a) (Subadditivity) For all states p and a, 

E{p®a) <E{p) + E{a). (4.16) 

(E5b) A special case of the above is when p and a are different numbers of the same 
state tensored together. For all states p, and m,n> 1, 

^(^®(m+n)-) < ^ E{p'^''). (4.17) 

(E5c) (Existence of a regularisation) For all bipartite states p, the limit 

E'^ip) = lim ' (4.18) 

n— >oo n 

exists. E°° is called the regularisation of E. This is the weakest of the additivity conditions, 
and is usually satisfied. 
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(E6a) (Convexity) Mixing of states does not increase entanglement. For all bipartite 
states p and a, 

E{\p + (1 - \)a) < \E{p) + (1 - \)E{a) (4.19) 

for < A < 1. 

Although (E6a) seems very natural as forgetting which of a set of prepared states you 
have should not increase entanglement, there are indications that the (asymptotic) distillable 
entanglement is not convex .SSTOlj . This is based on the fact that PPT bound entangled 
states together with NPPT bound entangled sates (if they exist) can be distilled. Thus, 
even a mixture of two states with zero distillable entanglement can be distilled. A weaker 
condition which is shown in |DHE,n2j to hold also for the distillable entanglement is that 
the convexity need only hold on decompositions into pure states. 

(E6b) For any bipartite state and any pure state realisation p = YliPi iV'j) {'4^i\-, Pi ^ 
and J2iPi = 1' 

Eip) < Y.p^Ei\iJi) ii^il). (4.20) 

i 

Like for pure states, asymptotic weak additivity is sufficient for weak additivity, so (E4a) 
and (E4b) are actually equivalent. The same goes for the subadditivity conditions (E5a) 
and (E5b). 

4.2.3 Ed and Ec as extreme measures 

In |DHE,n2j the exact conditions for Ed and Ec to be the lower and upper bound for an 
entanglement measure were investigated. Three versions of the theorem were found. With 
different conditions on the measures either the nonregularised or the regularised measure 
would be bounded by Ed and Ec- 

(i) For an entanglement measure E satisfying (Ela), (E2a), (E3a) and (E4b) (and thereby 
(E4a)), for all states p, 

Ed{p) < E{p) < Ec{p). (4.21) 

These conditions are very strong, and until recently no function was known that satisfied 
them. The "squashed entanglement" |CW04j . however, satisfies all conditions. 

(ii) For an entanglement measure E satisfying (Ela), (E2a), (E3a) and (E5c), then for 
all states p, the regularised version E^ (which exists by (E5c) and always satisfies (E4a), 
but not necessarily (E3a)) is bounded by Ed and Ec- 

Ed{p) < E^{p) < Ec{p) (4.22) 

These conditions are easier to satisfy, and they are satisfied by both the entanglement of 
formation and the relative entropy of entanglement. 

(iii) Let E be an entanglement measure satisfying (Ela), (E2a), (E3b) and (E6b). Then 
if weak additivity (E4a) holds, 

Ed{p) < E{p) < Ec{p) (4.23) 

and if subadditivity (E5a) holds, 



Ed{p)<E'^{p)<Ec{p)- 



(4.24) 
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4.2.4 Relative entropy of entanglement and other distance based mea- 
sures 

Along with the first attempts to give reahstic conditions for entanglement measures jVPRKOTl 
|VP98], aclass of measures that satisfied the conditions was introduced. The conditions con- 
sidered was that the measure would be zero if and only if the state is separable (EOa), that 
local unitary operations leave it constant (E2c) and that its expectation value does not 
increase under LOCC (E2a). 

The class of measures is based on some distance function (but not necessarily a metric) 
on the set of density matrices. The measure is then the distance from the state in question 
to the nearest separable state. That is, for a distance function D{p,a) an entanglement 
measure corresponding to the distance function could be defined as 

where S{Tl) is the set of separable states on 7i. The conditions imposed on E is now 
converted into conditions on D. 

The quantum relative entropy^ can be used as a distance function on the set of density 
operators. It is defined as 

Sip\\a) = tT[p{logp-loga)]. (4.26) 

The relative entropy is not a metric and is not even symmetric. It is nonnegative, and 
zero only for identical density operators. The same unitary operation on both states leaves 
it invariant. The quantum relative entropy can be interpreted as a distinguishability of 
quantum states. More precisely, suppose we are given a large (but finite) number n of 
quantum systems that are all in the same state, which is either p or a. Our task is to perform 
measurements to infer which. The probability for inferring from optimal measurements on 
the composite system, that the given state is p when it really is a is |VP.TK97; 

P„(a ^ p) = 2-"^WI'^) (4.27) 

(for large n). This expression is not symmetric in p and a. 
To illustrate this, consider the qubit states 

p=|0)(0| and a=^|0)(0| + i|l)(l|=^l2 

i.e. one arbitrary pure state (denoted |0)) and the completely mixed state. These states are 
optimally distinguished by measurement if we measure each particle in the {|0) , |1)} basis. 
In the case that we are performing measurements on a, the probability for measuring |0) and 
|1) are both |. If all measurements give |0), we would wrongly conclude that the state was 
p. The probability for that to happen is 2"". Conversely, if we perform the measurements 
on systems in state p, we will always get the result |0), so the probability for confusing it 
with a is zero. 

^For reviews of the role of the relative entropy in quantum information theory, see |Ved02lfSW02| . 



(4.25) 
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The quantum relative entropy can be calculated from the eigenvalues and eigenvectors of 
the density operators as shown in Appendix In the above example, it can be calculated 
from (|A.4() giving ^(pIIit) = 1 and ^(fTll/)) = +oo. This is consistent with the probabilities 
above and H4.27|) . It is a general feature that S'(cjIIp) = +oo when /? is a pure state. This is 
because for any pure state there is a complete measurement for which the outcome is 100% 
certain, and using this measurement would make it impossible to confuse it with another 
state. 

This property alone makes the relative entropy with the arguments reversed, D{p,a) = 
S{a\\p), unusable in (|4.25j) as a distance function. It would give infinity for any pure en- 
tangled p. S{p\\a) on the other hand, behaves well for separable a and the entanglement 
measure generated by it, 

Erip)^ mi S{p\\a), (4.28) 

creS{H) 

is called the relative entropy of entanglement |VP98j . It could be said that it quantifies 
the unlikelihood for a separable state to give measurement outcomes consistent with the 
entangled state. 

From its definition and the fact that ^(pIIit) = only when /? = o", it is obvious that 
the relative entropy of entanglement satisfies (EOa) {Er{p) = iS p G S{7i)). It reduces to 
the entropy of entanglement on pure states |VP98j , and therefore satisfies the normalisation 
criterion (Ela) and continuity on pure states (E3b) by the uniqueness theorem. It is nonin- 
creasmg under LOCC operations (E2a) |VPRK97| . Continuity (E3a) was shown in |DH99| . 
There have been found counterexamples to weak additivity (E4a) |VW01j and thereby also 
to asymptotic weak additivity (E4b). Actually, for very high dimension of a certain class 
of states, Er{p) ~ Er{p ® p). It can also be shown that it is both subadditive (E5a) and 
convex (E6a) |DHH,02| . 

The regularisation (E5c) 

Er(p) = E^ = hm ' (4.29) 

n— >oo n 

of the relative entropy of entanglement satisfies automatically (EOa), (Ela), (E2a), (E5a) 
and (E6a) as those conditions extends to any regularisation of a measure that satisfies them 
|DHR02j . In addition, as any regularisation does, it is satisfies weak additivity (E4a). 

Any distance based measure, including the relative entropy of entanglement, can be 
varied by considering other classes of states than the separable ones. What is important is 
that the class is closed under LOCC. For instance, the set can be extended to PPT states. 
This will give a measure which is easier to calculate and inferior or equal to the version for 
separable states. For 2x2 and 2x3 dimensional systems the measures coincide as the PPT 
is then equivalent to separability. 

Even though the quantum relative entropy is by far the most used of the distance based 
measures, other distance functions provide usable entanglement measures as well. Another 
measure of distinguishability of quantum states is the Bures metric (Bur69. ■Hiib92. .FC95j , 
in this context defined as |VP98j 



DB{p\\a) = 2-2^/F{p,a) 



(4.30) 
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where 

F{p,a)^[tv{{V^pV^y/'}]' (4.31) 

is called Uhlmann's transition probability |Uhl76j . The Bures metric is a true metric and 
thus symmetric. It is a generalisation of the pure state Fubini-Study metric to mixed 
states |Hub92j . The entanglement measure generated by it satisfies (EOa), (Elb) and (E2a). 
Unlike the relative entropy of entanglement, however, it doesn't reduce to the entropy of 
entanglement on pure states. Actually, it is smaller than the entropy of entanglement 
|VP98j. 

There may be other distance functions that generate entanglement measures with nice 
properties, but the two mentioned are the ones used in practice. The Hilbert- Schmidt 
distance 

Dhs{p, ct) ^ \\p - aWjjs = tr [(p - af] (4.32) 

was suggested as a candidate |VP98|IWT99j . However, it was shown by Ozawa lOzaOflj that 
the distance violated one of the criteria it was conjectured to have (nonexpansion under 
physical operations), so it did not satisfy the sufficient conditions for it to be nonincreasing 
under LOCC. On the other hand, the very similar trace norm distance 



DTip,a) = \\p-a\\i=tT [V{p-'J?\ (4.33) 

has been shown to generate an entanglement measure which is monotonic under LOCC 
|EAPn3| . 



4.2.5 Entanglement of formation 

The entanglement of formation was historically the first entanglement measure to appear 
|BDSW96] . In the literature it is occasionally called entanglement of creation. It was meant 
to be the asymptotic entanglement cost, but this interpretation rests on the weak additivity 
which is strongly conjectured, but still not proved. However, the regularisation of it has 
been shown to be equal to the asymptotic entanglement cost |HHTf)lj . 

The entanglement of formation is a straightforward generalisation of the entropy of 
entanglement to mixed states. Remember that in the asymptotic limit, the entanglement 
cost in Bell states of preparing a pure state is given by the entropy of entanglement. Thus, 
it was natural to define the entanglement of formation for a pure state as the entropy 
of entanglement. For a given ensemble of pure states £ = {pi, the entanglement of 
formation is the average of the entropy of entanglement for the states in the ensemble 

i^/(f)^ J^p.i^ijdV',)). (4.34) 

i 

A mixed state can be realised by a multitude of pure state ensembles, with different en- 
tanglement of formation. As any of those ensembles realises the mixed state, the natural 
definition for the entanglement of formation for a mixed state is the entanglement of forma- 
tion for the "most economic" ensemble. That is, the entanglement of formation for a mixed 
state is defined as 

Ef{p)^m^Y,P^EEm) (4.35) 

i 
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where the infimum is taken over all ensembles e = {pi, iV'j)} that realises the state p. 

How come entanglement of formation is not necessarily equal to the entanglement cost? 
It is easy to see that the asymptotic entanglement cost cannot be greater than the entan- 
glement of formation. Imagine that we are to produce a large number riout (f^out — > oo) of 
some bipartite state /Oout- We find the pure state realisation {pj, of po^t that has the 
lowest average entropy of entanglement (equal to the entanglement of formation) . Then we 
produce this ensemble by producing the pure states in amounts corresponding to the prob- 
ability distribution {pi}. In the asymptotic limit this can be done reversibly. When we mix 
the produced states (discard the information saying which is which) , we have a collection of 
systems in the state pont- Since the entropy of entanglement is conserved in the first step, 
and Ee = 1 for Bell states, we have produced systems in state /?out at a cost equal to the 
entanglement of formation. 

The question is if we can do better than this. Intuitively the above is the most economical 
way, as any other realisation of pout would give a higher or equal cost. But there is another 
way of producing a large number of pout- We can produce Hout/N copies of the state p^^- 
It is not known whether the optimal pure state realisation of this state will have an average 
entropy of entanglement times that of the state Pout- In other words it is not known if 
Ef is weakly additive, Ef{p) = ^Ef{p'^^). If it is, then Ej = Ep = E^ which is equal to 
the asymptotic entanglement cost Ec- Additivity has been shown for some special classes 
of states |HNm I lYWrnl IV I )( ^021 IShin2j . but whether it holds in general, is one of the big 
open questions in quantum information theory. 

From the definition of Ef^ it is easy to see that it satisfies (EOa) and (Ela). (EOa) by 
the fact that the optimal pure state decomposition of a separable states is into a mixture 
of product states, (Ela) by the fact that it reduces to the entropy of entanglement on 
pure states. LOCC monotonicity (E2a) was shown when it was introduced |BDSW96| . and 
continuity (E3a) was shown in INieODj . Subadditivity (E5a) follows from the discussion of 
additivity above, one may use the same type of decomposition for a tensor product, but one 
might be able to do better, which would give a lower Ef. Convexity (E6a) follows more or 
less directly from the definition. 

The way the entropy of entanglement was extended to mixed states in the case of the 
entanglement of formation, can also be used for other measures on pure states. Measures 
constructed this way are called convex roof measures |TJh]98llHorfrrj . The method extends 
a measure of some set (here pure states) to the convex hull^ (here mixed states), where it 
is the largest function that is convex and compatible with the measure on the original set. 
The method was used in |L(X)Kn3] construct an entanglement measure that extends the 
entanglement measure called negativity to mixed states in another way than the original 
measure. 

A concept related to the entanglement of formation is the concurrence |IIW97| IWoo98j . 
It is defined for a system of two qubits. For a general state p of two qubits, let p be the 
spin-flipped state 

p = {Y®Y)p*{Y®Y) (4.36) 

where the Y is the Pauli Y operator and p* is the complex conjugate of p, both taken in 

*The convex hull of a set is the set of all elements that can be written as a convex combination of the 
original set. For instance, the convex hull of two points is the line connecting them. 
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the standard basis {|00) , |01) , |10) , Let the Hermitian matrix R be defined as 



R^^JVppVp (4-37) 

with eigenvalues in decreasing order {Aj}. The concurrence is then defined as 

C{p) = max{0, Ai - A2 - As - A4}. (4.38) 

The concurrence is monotonic under LOCC, and can thus be used as an entanglement 
measure for two qubits. The great advantage is that it is easily computable. But more 
important is that it is directly related to the entanglement of formation, providing an explicit 
formula for the entanglement of formation in the case of two qubits. Let the function £ be 
given by 

£(C,..(i±4H^) (4.39) 

where h is the binary entropy function 

h{x) = — xlogx — (1 — x) log(l — x). (4.40) 

Then the entanglement of formation is simply given by |HW97l IWoo98j 

Ef{p)=£{C{p)). (4.41) 

There have been various attempts to generalise the concept of concurrence. Uhlmann 
|UhlOOj considered general conjugations instead of the special (|4.36|) . This was further 
generalised to concurrence vectors by Auderaert et al. |AVD01) . The previously mentioned 
convex roof extended negativity coincides with the concurrence in the case of two qubits 
|LCOK03] . and thus provides another generalisation. Some of those generalisations - along 
with other aspects of concurrence and entanglement of formation - are reviewed in |Woo01j . 



4.2.6 Negativity 

The entanglement measures discussed so far have all had a serious drawback. Their definition 
includes some kind of optimisation, which make their evaluation very difficult. In fact, 
expressions for the relative entropy of entanglement and entanglement of formation has 
only been calculated for highly symmetric states (e.g. |VP98I ITVOni IVWnil IAI)VWn2j ^ 
and low dimensional cases (e.g. |HW97l rWoo98j ) . For general states, evaluating a measure 
includes heavy numerical calculations. With the aim of introducing a computable measure 
of entanglement, two related quantities were defined by Vidal and Werner |VW02aj . They 
can both be seen as a quantification of the PPT criterion for separability, and because of 
that they have the disadvantage that they fail to distinguish between separable states and 
entangled PPT states. Both quantities are based on the trace norm of the partial transpose 
(in some basis) of a state, which can easily be calculated using standard linear 

algebra packages. The first quantity is the negativity 




(4.42) 
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This quantity is equal t o the abs olute value of the sum of negative eigenvalues of p^^ . It 
was first introduced in ZHSLQsl Appendix B] (without the factor i) and later shown to 



not increase under LOCC |VW02aj . The other quantity is the logarithmic negativity 

EM{p)^\og\\p^^\\u (4.43) 

which is not strictly monotonic under LOCC, but does not increase for a subclass of LOCC. 

For all separable states (and more generally all PPT states), both M and Ej^ vanish. 
This is easily seen from the fact that the partial transpose operation does not change the 
trace of the density operator. If p^^ is a positive operator (i.e p satisfies the PPT criterion 
for separability), then ||/o'^^||i = 1, since is the sum of the absolute values of the 

eigenvalues of p^^ and all eigenvalues are positive. For states violating the PPT criterion, 
IIp^^IIi > 1. In this sense, both quantities quantify how much the state violates the PPT 
criterion. 

A lot of properties of M and Ej^ were deduced in |VWn2aj . As already mentioned, 
N is zero for all separable states, but is also zero for PPT entangled states. Hence, it 
satisfies (EOb) but not (EOa). It is ^ for Bell states, so it doesn't satisfy the normalisation 
criteria either. It could be made to satisfy (Elb) by defining it as the double value, but the 
normalisation wouldn't extend to maximally entangled states in dimension d x d (Ela). It 
is monotonic under LOCC (E2a), convex (E6a), and superadditive, satisfying M{p(d cr) = 
M{p) + J\f{a) + 2J\f {p)J\f [a) . Thus, neither additivity of any kind (E4a) nor subadditivity 
(E5a) is satisfied. 

Like AA, logarithmic negativity Ej^ satisfies (EOb) and not (EOa). However, unlike N it 
also satisfies the normalisation condition (Ela). But it does not coincide with the entropy 
of entanglement for pure states in general. Actually Ej^ > Ee on pure states, with equality 
only for maximally entangled states. It is not monotonic under LOCC (E2a), but is an 
upper bound for the distillable entanglement. Nor is it convex (E6a). The most appealing 
property of Ej^f is probably that it satisfies (strong) additivity, Ej^f{p®a) = Ej\f{p) + Ej^{a). 
From this, weak additivity (E4a) and subadditivity (E5a) comes for free. 

The negativity can be generalised to give rise to other entanglement measures. This is 
done by replacing ||/0"^^||i by another norm ||p||. The generalisation comes naturally if we 
define the trace norm in another, but equivalent way. Any Hermitian matrix can be written 
as a difi'erence of two positive operators, 

A = a+p^ - a^p~ (4.44) 

where p^ are density matrices and a± are nonnegative numbers. There exists a decompo- 
sition of the form (|4.44l) which minimises a+ + a_ . For this decomposition the trace norm 
is ll^lli = a+ + a_, and a_ is the absolute sum of the negative eigenvalues of A |VW02aj . 
From this we get another expression for the negativity, namely 

M{A) = inf{a_ I A^B = a+p+ - a^p'}. (4.45) 

Minimising a_ is in this case the same as minimising a+ + a_ , since by taking the trace of 
the condition, we get o+ = o_ +tr(^^s). If j^Tb 

is a density matrix, A is PPT, and we can 
take a_ to be zero, so the negativity vanishes as it should for PPT states. 
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The condition A^^ = a+p"*" — a_p~ restricts the set of density matrices over which 
we can take the infimum. Another norm and another negativity can be obtained by taking 
another set S from which must be picked. Norms that can be defined in this way 
corresponding to a compact set S are called base norms. The norm corresponding to a 
general compact set S is then 

\\A\\s = inf{a+ + a_ \ A = a+p~^ - a-p~, a± >0,p^e S} (4.46) 

and the corresponding negativity 

J\fs{A) = inf{a_ I A = a+p+ - a^p-, a± > 0, G S}. (4.47) 

We get back the usual trace norm and negativity from this by taking S as the set of all 
Hermitian (not necessarily positive) matrices with unit trace such that the partial transpose 
is positive. In mathematical terms, S = {A \ A = A^ , tr(^) = 1, A'^^ > 0}. 

By choosing S as the set of all separable states, another entanglement measure arises. 
It was described prior to the introduction of negativity and was introduced in |VT99j . It 
is called robustness of entanglement and was defined as the minimum amount of mixing 
with locally prepared states needed to make the state separable. Also other measures have 
been related to negativities |Horf)lj , such as the measures introduced by Rudolph based on 
the greatest cross norm |Rud01j . It should be noted that only the original negativity has 
the advantage of being easily computable. The other negativities need to be evaluated by 
finding the infimum over a large parameter space. 



4.2.7 Squashed entanglement 

Of the entanglement measures considered so far, only few are additive, and when they 
are some of the other greatly desired conditions are not fulfilled. Quite recently a new 
entanglement measure was proposed, having more nice properties than its predecessors. It 
is motivated from the so called intrinsic information of classical cryptography and is called 
squashed entanglement |CW04j . It is defined as 

^sq(p) = inf |^J(A; B\E) \ extension of p to We| (4.48) 

The infimum is taken over all extensions to a third subsystem E, so p = p^^ = tiE{p^^^)- 
I(A;B\E) is the quantum mutual conditional information defined as 

I{A; B\E) = S{p^^) + S{p^^) - S{p^^^) - S{p^). (4.49) 

Most of the properties of Esq were shown along with its introduction as an entanglement 
measure in |CW04j . The squashed entanglement vanishes on all separable states (EOb), but 
it is not known whether it also vanishes for some entangled states (EOa). If there exists 
an extension to 7{e with dimTiE < oo for which I{A;B\E) vanishes, it implies that p^^ 
is separable. But the infimum in (|4.48|) is taken over all extensions, and thus Esq{p^^) 
may be zero even if no finite dimensional extension makes I {A; B\E) zero. It coincides with 
the entropy of entanglement on pure states, and thus satisfies the normalisation condition 
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(Ela). It is also bounded from below by the distillable entanglement and from above by the 
entanglement cost. As a proper entanglement measure it does not increase under LOCC 
(E2a) and it is convex (E6a). When it was introduced, it was shown to be continuous on 
most of the state space, and indications that it was continuous on the whole state space were 
presented. Continuity was finally proved by Alicki and Fannes |AF04j (E3a). The measure 
is strongly additive, making it satisfy both weak additivity (E4a) and subadditivity (E5a). 
So except from (EOa), which is still unknown, squashed entanglement satisfies the strictest 
versions of all conditions presented. 

In |VW01j various additivity properties were discussed. The additivities we have pre- 
sented so far are for independent pairs, described by tensor products. But if there are 
entanglement or classical correlations between the pairs as well, they may be taken advan- 
tage of. Consider the state p'^^^ ^ which lives on the Hilbert space Ha ^Hb® Ha' 'SiH-b', 
where Ti.A ® Ti-A' is controlled by Alice and TIb ® Hb' by Bob. Each of the primed and 
unprimed subsystems are individually described by tracing out the other, 

= tr^,^,(/^^'^'), p^'^' = tr^^(/^^'^'), (4.50) 

and if they are independent p"^^^ ^ = p^^ ® p"^ ^ . For a strongly additive entanglement 
measure this implies E{p^^^ ^ ) = E{p^^)+E{p^ ^ ). But if there are correlations between 
the pairs that can be exploited, we would expect 

E^pABA'B'^ > Eip^"") + (4.51) 

This property is called strong superadditivity |VW01j and is satisfied by the squashed en- 
tanglement |('Wfl4j . 

Like most other entanglement measures, the definition of squashed entanglement involves 
an optimalisation which makes the measure hard to compute. Due to its recent introduction, 
not much work has been published about its calculation on special sets of states, nor about 
efficient numerical algorithms. One of the important open questions is whether it vanishes 
on any entangled states, especially on the bound entangled ones. 

4.2.8 Other entanglement measures 

In addition to the entanglement measures discussed, other measures have been proposed 
without g aining the same importance. One of them is the entanglement of assistance |Coh98l 
IDFM+ 99] . It is based on the following idea. Many copies of an entangled tripartite pure 
state, l^'"^'^*-^), is shared between Alice, Bob and Charlie. Now, Alice and Bob want to 
use the state to perform some task involving the use of bipartite entanglement. If they 
ignore Charlie completely, they share the state p^^ = trc(|^'^^*-^) (^'^^'-^|). In other words 
they share a mixed state, and it might even be separable. But Charlie wants to help, 
and Alice and Bob may communicate classically with him. Therefore, Charlie can perform 
local operations on his part of the system and tell Alice and Bob about the measurement 
outcomes. A measurement by Charlie will then leave the subsystem of Alice and Bob in some 
pure state, which depending on the measurement outcome may or may not be entangled. 
Alice and Bob can then discard the pairs that are not entangled, and turn the rest reversibly 
(in the asymptotic limit) into Bell states. If each measurement by Charlie leaves them with a 
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state l^^i) with probability , the number of Bell states produced per copy of the tripartite 
state will be equal to the average entropy of entanglement of the ensemble £ = {pi, \ 
namely 

E{£) =Y,P^EE{\'i'f''))■ (4.52) 

i 

Now, Charlie wants to do the measurements in a way that leaves Alice and Bob with as 
much entanglement as possible on average. It turns out that the optimisation can be taken 
over all ensembles £ = {pi,\"^f^)} consistent with p^^. Thus, what Charlie can achieve 
does not depend on the initial tripartite state \^^^'^), only on p^^ , so all pure tripartite 
states for which p^^ is the same, allow Charlie to help Alice and Bob to the same result. 
The entanglement of assistance is the average entropy of entanglement Alice and Bob are 
left with after Charlie has applied his optimal strategy, namely 

Eaip^^) = supY.p.EEil^f'')). (4.53) 

^ i 

We can see that this is dual to the entanglement of formation (|4.35j) , which is the same with 
an infimum instead of the supremum. Even though the entanglement of formation is believed 
to be additive, there are states for which the entanglement of assistance is superadditive, 

£^a(/^ ^^^) > K(/^) + Eaia""^). (4.54) 

A recently introduced measure related to the entanglement of assistance is the localizable 
entanglement |VPC04j . It is really a multipartite measure for a large number of qubits, but 
the goal is, just as for the entanglement of assistance, to optimise (localise) the entanglement 
between two subsystems by performing measurements on the others. If we had allowed global 
measurements on the other qubits, we could have regarded it as a single "Charlie" system, 
and it would have reduced to the entanglement of assistance. Because we can only perform 
local measurements on each individual qubit, the localizable entanglement is necessarily 
smaller than the entanglement of assistance. 

Another entanglement measure recently introduced is the witnessed entropy of entangle- 
ment |BV04j . It is based on the concept of entanglement witnesses |HHH96j . An entangle- 
ment witness for an entangled state o" is a Hermitian operator W for which tr{Wa) < and 
tr{Wp) > for all separable states p. We conventionally normalise entanglement witnesses 
to have tr(VF) = 1. We say that an entanglement witness is optimal for a state a if 

tr{Waa) < tT{Wa) (4.55) 

for all entanglement witnesses W. The definition of the witnessed entropy of entanglement 
is based on the optimal witness, 

EM = log{D/d) max[0, - tT{Wpp)] (4.56) 

where D is the dimension of the total Hilbert space and d the smallest of the dimensions 
of the Hilbert spaces of the subsystems. Note that the definition applies equally well to 
bipartite and multipartite systems. It is not additive, but satisfies the rest of our condi- 
tions for entanglement measures, (EOa), (Ela), (E2a), (E3a), (E5a) and (E6a). The main 
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advantage of this measure is that it can be approximately calculated for all mixed states. 
Finding the exact optimal entanglement witness, is a hard optimisation problem, but using 
approximation techniques, there are algorithms to find the E^j to a given precision. 

Other measures can be defined operationally in the following way. Take your favourite 
application of entanglement, study how well or bad it performs for diff'erent sets of states, 
define an entanglement measure from this and check that it satisfies at least some of the 
conditions for entanglement, above all LOCC monotonicity. One of the operational measures 
that have had limited impact on the community is an entanglement measure introduced 
by Hiroshima |Hir02j based on the capacity of dense coding |BW92j . Another operational 
measure by Biham et al. |BN()n2j is derived from Grover's search algorithm |(-rro96|r(Tro97j . 
It is based on the fact that the algorithm performs worse the more the input state is 
entangled. 



Chapter 5 



Beyond bipartite entanglement and 
finite dimension 



The systems we have considered until now have all been bipartite systems with finite dimen- 
sional Hilbert spaces. Some of the concepts and measures considered are valid also in the 
multipartite systems and for entanglement of continuous variables. Others lose their mean- 
ing in this more general setting. In this chapter we consider some of the steps that need to 
be taken to generalise the concepts from the previous sections to multipartite systems and 
systems with continuous variables. 

5.1 Multipartite systems 

While the basic structure of bipartite entanglement is well understood, multipartite entan- 
glement is still an active field of research. In bipartite systems the pure state entanglement 
in the asymptotic regime is well understood, since the von Neumann entropy of the reduced 
density matrices can be conserved under transformations. In general the Schmidt decompo- 
sition cannot be extended to multipartite systems |Per95j . so the eigenvalues of the reduced 
density matrices of each system need not be the same, nor do the von Neumann entropies. 

The questions about entanglement literally acquire new dimensions when more that two 
subsystems are considered. Between two parties all relations are between those two parties, 
and in terms of quantum communication there is only one single channel. Once we increase 
the number of parties to three, we have three pair relations in the system. Still, those three 
pair relations do not describe all correlations in the system (neither classical nor quantum). 
Consider for instance the pure tripartite entangled state called the GHZ^ state 



If we trace out any one of the systems, we get the density matrix shared by the two remaining 
parties 



GHZ) = ^(|000) + |111)). 






^Daniel M. Greenberger, Michael A. Horne and Anton Zeilinger 
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which is separable. Thus, the state contains no bipartite entanglement between the systems, 
but the tripartite entanglement is maximal. 

5.1.1 Separability 

Some notions can be taken directly over from bipartite systems, though. The generalisation 
of separability from ()3.3|) is straightforward. A state is separable if and only if it can be 
written as a convex combination of product states; 



with Pi > and X^jPi = 1. 

The classification of states does not end with this. An entangled multipartite state may 
or may not still be entangled after one of the subsystems are ignored (i.e. traced out). A 
multipartite state that loses its entanglement when any of the subsystems is traced out, is 
called multiseparable |Tha99j . 

The PPT criterion for separability trivially extends to the multipartite case. If a multi- 
partite state is separable, then transposing any number of the subsystems must also give a 
valid (i.e. positive) density operator. 

5.1.2 Measures of entanglement 

The definition of entanglement cost and the distillable entanglement is based on a standard 
bipartite state, namely the Bell state. This is a bipartite state, so it is not obvious how to 
extend these measures to the multipartite case. Nielsen |Nie01j argued that this could by 
overcome by generalising the definition from Bell states to a definition in terms of qubits 
transferred (in the case of creation) or the number of CNOT gates that could be implemented 
(in the case of distillation). For instance, to create a Bell state it is necessary to transfer one 
qubit. Likewise the transfer of one qubit can be simulated by teleporting a state, spending 
one Bell state. If the entanglement cost is defined in terms of how many qubits that need 
to be transferred, it will coincide with the usual definition for bipartite states. 

Most of the conditions we imposed on measures in section 14.2.21 generalise to the multi- 
partite case. The only possible exceptions are of course the normalisation conditions (Ela) 
and (Elb) which are based on Bell states. Those would be unnatural to impose on mul- 
tipartite entanglement, but it should be possible for a measure to satisfy them for states 
where only two of the subsystems are entangled. 

The relative entropy of entanglement and other distance based measures do not have 
any reference to the bipartite case in their definition. Therefore, they are just as valid in 
the multipartite case as in the bipartite case. 

The entanglement of formation suffers the same problem as the entanglement cost, and 
can be generalised in the same way. For a pure state we define the entanglement of formation 
as the number of qubits we need to transfer in order to prepare it (LOCC comes for free). 
Then we extend the definition to the mixed states by the convex roof method as usual. 
Another more formal generalisation was suggested by Wang jWanOlj . 




(5.3) 
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The meaning of entanglement measures in the multipartite setting is more vague than 
in the bipartite. In the bipartite setting it quantifies in some way how much that can be 
achieved by the two parties together. In the multipartite case there are more possibilities. 
For instance there may be some entanglement between all of the parties, as in the state 

|VF)^-^(|100) + 1010) + 1001)) (5.4) 

which on tracing out one of the subsystems gives the density matrix 

p=^|00)(00| + ^|vI/+)(^+| (5.5) 

which is an entangled state with negativity M = 0.206 and entanglement of formation 
Ef = 0.550 (calculated from the expressions in sections 14.2.61 and 14.2. 5() . Or there may be 
no entanglement between the pairs, as in the GHZ state mentioned earlier. Both states can 
be prepared by sending two qubits (after one of the parties has created the state locally). 
While the different bipartite measures place different orderings on the states |VP00j , ordering 
is even more ambiguous for multipartite states. The types of entanglement are qualitatively 
more different, and a state being very much entangled in one sense (as defined by one 
measure), can have very little or no entanglement in another sense. All this says is that, 
since a multipartite system is a more complicated structure, we cannot expect one single 
parameter to be as good a description of the entanglement in it as it is in the bipartite case. 



5.2 Continuous variables 

Most of quantum information theory has been formulated in the context of finite dimen- 
sional Hilbert spaces. However, from an experimental point of view, entanglement is just 
as interesting in the continuous variable regime. This has led to efforts to quantify en- 
tanglement also in the case where the Hilbert spaces have infinite dimension. Two papers 
by Eisert et al. |ESP02| IEP03j summarise many of the new problems that arise and their 
possible solutions. 

First of all, when the dimension of the Hilbert space is infinite, we have states with 
an infinite amount of entanglement. This should come as no surprise, as the entropy of 
entanglement for a maximally entangled state m. dxd dimensions is logc? (c.f. criteria (Pla) 
and (Ela)). What is more worrying is that the set of pure states with infinite entropy of 
entanglement is actually dense in the trace norm on the set of pure states |ESP02j . This 
means that arbitrarily close to any product state is a state which has infinite entanglement. 
The solution to this problem is to consider only states where the mean energy is bounded 
from above. This is a reasonable physical assumption, and on this subset of states the 
entropy of entanglement is continuous. Other measures can also be defined on this subset 
without exhibiting strange behaviour. 

The definition of separability in a continuous variable bipartite system must be slightly 
altered from the one in section TS. 2. 11 The expression is still the same as (|3.3|) . namely 



i 



(5.6) 
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with Pi is a probabiUty distribution. But in infinite dimension we can only require that the 
state can be approximated arbitrarily well by an expression of the form ()5.6() . 

The PPT criterion for separability is also valid in infinite dimension and for the special 
class of so-called Gaussian states, it is also sufficient for separability |SimOOj . 

It is a well known fact in the finite dimensional case that there is a neighbourhood around 



the completely mixed state I/d where all states are separable |iZHSL98, PBOO^ Myr04 
However, the volume of the separable ball diminishes as the dimension of the system grows. 
In the infinite dimensional setting, even when the energy is bounded from above, it has been 
shown |ESP02j that arbitrarily close to any state (in the trace norm), there is an entangled 
state. 

The class of so-called Gaussian states has a special role in entanglement theory of con- 
tinuous variables. These are states with a Gaussian Wigner function (a quasi-probability 
distribution in phase space which completely describes a quantum state). These states are 
often encountered in the laboratory, and they can be completely described by two param- 
eters. For Gaussian states, and operators that takes Gaussian states to Gaussian states, 
many of the notions from finite dimensional entanglement theory have similar counterparts. 
See e.g. (EPOHl IK(;L( lGWK+n3l IWGK+n4l iREfUj. 



Appendix A 



Expressions for the quantum 
relative entropy 

The quantum relative entropy is defined as 

S{p\\a) = tT{plogp- ploga}. (A.l) 

We want to calculate it in terms of the eigenvectors and eigenvalues of p and a. We denote 
the eigenvalues of p (a) by Vi (si) and its eigenvectors by \ri) {\si)). In the case of degenerate 
eigenvalues we assume that the eigenvectors have been properly orthogonalised. 
The linearity of the trace gives 

S{p\\a) = tr{plogp} - tr{plogcr} 

where the first term can be identified as —S{p), where S{-) is the von Neumann entropy. 
This depends only on the eigenvalues of p: 

tr{/9 log p} = ^ rj log ri. 

i 

The second term is more complicated as the eigenvectors of the density operators in 
general are not the same. We decompose the density operators using the spectral decom- 
position. 

p log (7 = ( ^ ri In) {n I ) log ( ^ Sj \sj) {sj I ) 

i j 
i j 

= ^nlog{sj){ri\sj) In) {sj\ 

Then we decompose \ri) in the eigendirections of a, \ri) = (sfc|rj) \sk)- 
ploga = ^rilog(sj) {nlsj) {skln) \sk) {sj\ 

ijk 
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Now we can take the trace of this to get the second term 

tr{/>logcr} = ^rjlogSj \ {ri\sj)f 
ij 

We can now put the two terms together to get the expression we wanted; 



S{p\\a) = J^r.j logr, -^^logsj \{n\s,)\^ |. (A.2) 



This may also be written as the single sum 



S{p\W) = ^ I rilogn - {si\p\si) log Si \. (A.3) 

i 

In the case when the eigenvectors for the two density matrices are the same (i.e. p and 
a commute), we can arrange the eigenvectors so that {ri\sj) = 6ij. ()A.2|) then reduces to 

S{p\\a) = Y,ri{logn- log Si). (A.4) 
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